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Abstract. The zodiacal foreground for a highly sensitive space infrared interferom- 
^ : eter is predicted for various observing locations. For the predictions we use a model 

that was derived from measurements of the Cosmic Background Explorer (COBE). 
We find that at a wavelength of 10 fim 96% of the sky is darker than 1 AIJysT~^ for 
observations in the ecliptic plane at 5 AU, and 83% is darker than 0.1 MJysr~^. At 
1 AU, however, always more than 50% of the sky are brighter than IMJysr"^, even 
if the observations are made from 30° or 60° of latitude above the ecliptic plane, at 
10 or 20 ^m. Thus, according to the employed model, the foreground reduction by 
increasing the heliocentric distance of the observing location is more effective than 
\ by increasing the latitude. 
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1. Introduction 



o 

:^ 

2 ' The search for extra-solar planets (exo-planets hereafter for brevity) 

"t^ , and possibly for primitive forms of life on them has received much 

I attention since the mid-1990s and is pursued by NASA and ESA. In its 

studies ESA follows a proposal by Leger (1996) of a space infrared inter- 
^ I ferometer, for which initially five 1 m-class telescopes were envisioned. 

^ • More recent designs of the interferometer, which was named DARWIN, 

use six 1.5 m telescopes. At NASA a similar instrument with the name 
of Terrestrial Planet Finder (TPF) is under development. In addition 
to the search for exo-planets such an interferometer could also be used 
for general purpose astronomical imaging and spectroscopy. While this 
paper is intended to support the mission analysis for DARWIN and 
TPF, it is independent from the actual baseline mission scenarios. 

One of the main problems for the detection and analysis of Earth- 
sized exo-planets is the considerable amount of zodiacal infrared fore- 
ground radiation emitted by the interplanetary dust cloud. In the vicin- 
ity of the Earth either large telescopes have to be used (Angel, 1989) or 
the observations have to be integrated over long periods of time in order 
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to sufficiently suppress tfie zodiacal foreground contribution (ESA-SCI, 
200x). The alternative proposal by Leger is to place the interferometer 
at a larger heliocentric distance, where the zodiacal foreground in the 
lO-to-20 /nm range is reduced due to the lower interplanetary dust 
density and lower dust temperatures. There the zodiacal foreground 
becomes comparable to other sources of noise, like the radiation leakage 
from the central star. As an alternative to larger heliocentric distances 
it is conceivable that the zodiacal infrared radiation is reduced at higher 
ecliptic latitudes. In order to cover the proposed orbit options, we inves- 
tigate heliocentric orbits with aphclia at 1, 3, and 5AU and inclinations 
of 30° and 60° (Jehn and Hechler, 1997). 

In this paper we determine the zodiacal infrared emission at various 
locations (heliocentric distances and ecliptic latitudes) in the solar sys- 
tem using a model by Kclsall ct al. (1998). This model was fit to the 
data taken by the COBE infrared satellite. We will discuss our choice 
of the dust model below. 

The result of our analysis is the fraction of the sky where the zodiacal 
foreground is of the same order as the other noise contributors, which is 
NzL ~ 2 X lO'' photons hour~^. At 10 //m this translates to a maximum 
foreground surface brightness of ~ O.lMJysr"^ for 1.5m telescopes, an 
interferometric transmission of 20 %, a spectral resolution of 20, and a 
field of view opening angle of 3.4 arcsec. 



2. Zodiacal Dust Model 

Up to date, no infrared telescope has been flown beyond 1 AU or out of 
the ecliptic plane. We therefore have to rely on models to determine the 
zodiacal infrared foreground at the suggested locations of the interfer- 
ometer. If we assume that the dust emits blackbody radiation, the dust 
model has to provide the spatial distribution of the cross section and 
the temperature of the dust grains in order to determine the infrared 
brightness. 

There exist a number of interplanetary dust models that are derived 
from zodiacal light and infrared observations, as well as in situ mea- 
surements. Our approach here is to use the phcnomenological model by 
Kelsall et al. (1998), that was derived from the measurements of the 
COBE/DIRBE instrument. The COBE/DIRBE measurements repre- 
sent the most accurate sky survey at infrared wavelengths between 
1.25 jiia and 240 jiia so far. The Kelsall et al. model describes the 
interplanetary dust cloud as made up of three components: a smooth 
lenticular cloud, dust bands generated by collisions within asteroid fam- 
ilies (Dermott et al., 1994), and the Earth-resonant dust ring (Reach 
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et al., 1995). We only consider the dominant contribution from the 

smooth cloud. 

The drawback of our approach is that, since COBE was an Earth 
orbiting satellite, the model is only well constrained for predictions 
at 1 AU in the ecliptic. Thus the prediction of an observation at, 
for example, z = 1 AU above the ecliptic plane along a line of sight 
towards the ecliptic north pole is uncertain, because it is not known 
how much of the brightness measured by COBE at 2; = along the 
same line of sight is generated at AU < z < 1 AU, and how much at 
z > 1 AU. However, the COBE observations along lines of sight with 
different ecliptic latitudes do constrain the vertical profile given by the 
model, whereas the radial profile is constrained by lines of sight with 
different ecliptic longitudes. The use of the Kelsall et al. model has the 
advantage that wc don't have to make assumptions about the grain size 
distribution, because it is based on infrared observations which depend 
only on the spatial cross section density, not the spatial particle number 
density. The model also does not consider different grain temperatures 
for different grain sizes, it applies an average grain temperature that 
depends only on the distance from the Sun. 

We summarize the Kelsall et al. model. The Infrared brightness Ix 
along a line of sight (LOS) is given by 



h = [ dsn,Ms))Bx{T) (1) 
JhOS 
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where h is the Planck constant, c the speed of light, and k the Boltz- 
mann constant. The parameters To = 286K and 5 = 0.467 were fitted 
to the COBE data. Considering the offset and the tilt of the zodiacal 
cloud with respect to the ecliptic plane, the cross section density rics is 
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for 
for 



C < /U 



(3) 



z" = x' sin O sin z — y' cos f2 sin z + z' cos i 

The parameters are: ro = 1 AU, no = 1.13 x 10"''' AU~\ a = 1.34, 
/5 = 4.14, 7 = 0.942, fi = 0.189, i = 2.03°, n = 77.7°, xq = 0.0119 AU, 
yo = 0.00548 AU, and zq = -0.00215 AU. 

Note that we have not apphed a color correction and assume an 
albedo of and an infrared emissivity of 1. Figure 2 shows the infrared 
radiation emitted by the smooth interplanetary dust cloud, predicted 
by the Kelsall et al. model. At the shorter wavelengths the radiation 



power density [W AU"^ sr'' Hz '] 
<10-'' 10"' 10"' 10' 10' 10=>10' 



10 iJ.m / 




X [AU] 



Figure 1. Power distribution radiated by the interplanetary dust cloud according 
to the Kelsall et al. (1998) model at 10, 20, 60, and 100 /xm. The 2-direction is 
perpendicular to the ecliptic plane and the x-axis points towards the vernal equinox. 

is emitted from the region closer to the sun, whereas the more distant 
regions contribute more to the radiation at longer wavelengths. Thus, 
at 10 fj,m the zodiacal foreground is more localized. 
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3. Sky Visibility Prediction as Function of the Selected 

Orbit 

We predict the infrared brightness received from the zodiacal fore- 
ground along lines of sight distributed over the whole sky. To represent 
the direction of the line of sight we use two angles : the solar-relative 
ecliptic longitude Aecl — Aecl,© and the ecliptic latitude /3ecl- Figure 
3 shows sky maps of the zodiacal foreground brightness /j/ for in-ecliptic 
orbits at heliocentric distances of 1, 3, and 5 AU. Due to the localized 




Figure 2. Sky maps of the infrared surface brightness Iv of the zodiacal foreground 

at lOfim (a), (b), (c), and 20 /im (d), (e), (f). Panels (a) and (d) show the brightness 
at an in-ecliptic observing location at 1 AU, in (b) and (e) the observation is made at 
a heliocentric distance of 3 AU, and panels (c) and (f) show the brightness at 5 AU. 
The contour lines show limiting foreground brightnesses of 0.1 and 1 MJy sr"'^. 

nature of the 10 fim radiation, the installation of the interferometer 
at larger heliocentric distances is very effective in the reduction of the 
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foreground brightness at 10 /xm. The reduction is at least two orders 
of magnitude in most regions of the sky. The 20 ^m brightness is less 
affected, but still significantly reduced at larger distances. 

We show sky maps of the foreground at observing locations 1 AU 
from the Sun, at ecliptic latitudes of 30° and 60° in figure 3. According 




Figure 3. Sky maps of the infrared surface brightness I,j of the zodiacal foreground 
at 10 /um (a), (b), and 20 /xm (c), (d). Panels (a) and (c) show the brightness at an 
observing location at 30° ecliptic latitude and 1 AU, and panels (b) and (d) show 
the brightness at 60° ecliptic latitude. The contour lines show limiting foreground 
brightnesses of 0.1 and 1 MJy sr~^. 

to the Kclsall et al. model, the zodiacal infrared foreground for lines of 
sight towards the ecliptic north pole is reduced by one to two orders 
of magnitude, if an observing location 60° above the ecliptic plane is 
selected. The reduction factor is about the same at 10 and 20/um. Table 
3 shows the fraction of solid angle of the sky that is darker than 0.1 or 
1 MJy sr-^ at 10 and 20 /zm. 

4. Discussion 

We have predicted the zodiacal foreground brightness at observing loca- 
tions in the ecliptic at 1, 3, and 5AU, as well as at 30° and 60° above the 
ecliptic at 1 AU using the Kelsall et al. (1998) model. Since the model 
is constrained only by observations at 1 AU in the ecliptic plane, the 
predictions have to be taken with caution, especially at high latitudes. 
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Table I. Fraction of the sky with a zodiacal foreground brightness below 0.1 and 
IMJysr-^ 



fraction of dark sky fraction of dark sky 
heliocentric ecliptic (< 0.1 MJy sr"^) (<lMJysr"^) 



distance 


latitude'^ 


10 iJ,m 


20 /xm 


10 iim 


20 urn 


1 AU 


0° 


0% 


0% 


0% 


0% 


1 AU 


30° 


0% 


0% 


0% 


0% 


1 AU 


60° 


0% 


0% 


38% 


6% 


3 AU 


0° 


0% 


0% 


83% 


0% 


5 AU 


0° 


83% 


0% 


96% 


70% 



^ of observer's location 



The validity of the model within the ecliptic plane at larger heliocentric 
distances is supported by zodiacal light observations by imaging photo- 
polarimeter on board the Pioneer 10 spacecraft (Manner et al., 1976). 
From the Pioneer measurements it was concluded that the dust density 
decreases as r~'^ with heliocentric distance r, and the best-fit value for 
1/ is 1.0 to 1.5. This is in accord with the radial exponent of —1.34 used 
in the Kelsall et al. model. 

As an alternative to the Kelsall et al. model we could have used a 
model derived from in situ measurements, like the Divine (1993) model. 
The Divine model was fit to, amongst others, in situ data that was 
collected at high ecliptic latitudes by the Ulysses spacecraft. It has to 
be considered, however, that the in situ measurements are dominated 
by particles that are much smaller (0.1 to 1 /im) than those that dom- 
inantly contribute to the zodiacal infrared emission (10 to 100 /Ltm). 
Additional constraints to the high latitude population of interplane- 
tary dust in the Divine model come from radar meteor measurements, 
which allow the determination of meteoroid orbits. Radar meteors give 
valuable information about the inclination distribution of interplane- 
tary meteoroids. However, all observed radar meteors are caused by 
Earth-crossing meteoroids, which introduces a bias to the dataset. 

In summary it seems worthwhile to do the same prediction, that 
we have made with the Kelsall et al. model, with the Divine model 
for comparison. This would improve our confidence in the projected 
zodiacal brightness at the high-latitude observing position. We found, 
however, that the Divine model predicts an infrared brightness at 1 AU 
in the ecliptic that are about a factor of 4 too low compared with 
the COBE measurements along the same line of sight. Since the Divine 
model was also fit to observations by the Infrared Astronomical Satellite 
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(IRAS) (Staubach et al., 1993), and the brightness measured by IRAS is 
in agreement with the COBE measurements, wc conclude tliat an error 
in the fit procedure lead to wrong grain temperatures. Staubach et al. 
(1993) report temperatures of 187 K for macroscopic (100 /Ltm) grains 
at 1 AU, which is in deed very low. Only extremely reflective material 
would reach such a low equilibrium temperature at 1 AU. Kelsall et al. 
(1998) use an average grain temperature at 1 AU of 286 K. Whatever 
the solution to the problem is, the discrepancy of the values shown 
by Staubach et al. (1993) in their figures 3 and 4, and the numbers 
calculated with their equations (1) - (3) has to be resolved before the 
Divine model can be used to predict the infrared brightness. 



5. Conclusion 

Our calculations have shown that an orbit with an aphelion at around 
5AU offers a better sky coverage (i.e. percentage of the sky for which the 
zodiacal foreground is below a certain level) than an inclined orbit at 
1 AU. Even on orbits with inclinations of 60°, the infrared emission 
of the zodiacal light prohibits the detection of exo-planets in large 
parts of the sky. In order to reach such orbital inclinations enormous 
changes in the orbital velocity are required that are not feasible with 
conventional chemical propulsion systems and that are very costly with 
electric propulsion engines (Jehn and Hechler, 1997). Therefore, for in- 
frared interferometer missions that are affected by the zodiacal infrared 
foreground, like the proposed DARWIN mission, orbits in the ecliptic 
plane with an aphelion around 5 AU are more promising. 

The strong model dependence of our prediction for observing loca- 
tions outside 1 AU shows, that precursor missions dedicated to make 
infrared measurements between 1 and 5 AU are desirable. 

Finally it has to be mentioned that during this investigation the 
baseline for the DARWIN mission has changed. Now the interferometer 
is planned to be installed in one of the quasi-stable Sun-Earth libration 
points (probably L2), which lie 1.5 x 10^ km from the Earth on the Sun- 
Earth line, in order to increase the reliability and to simplify operations 
and communications (ESA-SCI, 200x). At this observing location the 
photon noise from the zodiacal foreground dominates all other sources 
of noise like leaking starlight, detector noise, and galactic background. 
As can be seen in figures 3 (a) and (d) we expect the zodiacal foreground 
to be in the order of lOMJysr"^ at 10 //m. This means that in order to 
reduce the signal to noise ratio to a level where a terrestrial exo-planets 
can be detected, the observation of the target star has to be integrated 
over a long period of time. 
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